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Abstract. The aim of this work is an analytic investigation of differential equations 
producing mirror maps as well as giving new examples of mirror maps; one of these 
examples is related to (rational approximations to) We also indicate certain 

observations that might become a subject of further research. 



The existence of this paper is due to the foUowing observations. With Apery's 
proof of the irrationahty of ^(2) and C(3) followed certain 2nd- and 3rd-order linear 
differential equations (see [Bel]-[Be3], [BP]). Doing a similar construction for ^(4) 
resulted in a 5th-order differential equation. This equation was similar to the linear 
differential equations occuring in Calabi-Yau theory (except the order was 5 instead 
of 4). Computing the Lambert series of an analogue of the Yukawa coupling we got 
integer coefficients divisible by the square of the degree (not by the cube as in the 
Calabi-Yau case). Then we managed to pull back the 5th-order differential equation 
to one of order 4, which had all the properties of a Calabi-Yau equation. This part 
is the main objective of Sections 1-4. Quite recently we found two more examples 
relating fourth order differential equations to simultaneous rational approximations 
to di- and trilogarithms (see the end of Section 4). 

Collecting known cases of the Calabi-Yau equations, 14 of those were classical 
hypergeometric and 15 remaining were discovered by V. V. Batyrev, D. van Straten 
et al in [BSl], [BS2], [Str], we tried to extend the list by application of the algorithm of 
creative telescoping due to Gosper and Zeilberger. This seems to be a pure machinery 
(cf. Section 5 below), in which we cannot predict any end. We also tried to figure 
out possible algebraic transformations between them and our new example related 
to C(4)- We did not succeed in this, but doing that we found two other Calabi-Yau 
equations emanating from quadratic transformations of hypergeometric Calabi-Yau 
solutions of type 4F3. This is the subject of Section 6. Section 7 describes another 
(rather rich!) possibilities for deriving fake 4th-order linear differential equations. 
Summarizing methods of Sections 6 and 7 resulted in new classes of transformations 
which we present in Section 8. 



■i-The work of the second author is supported by an Alexander von Humboldt research fellowship 
and partially supported by grant no. 03-01-00359 of the Russian Foundation for Basic Research. 
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The divisibility properties of the coefficients of the Lambert series for the Yukawa 
couphng K{q) are equivalent to high-order Kummer congruences between the coef- 
ficients of the power series for K{q). The latter coefficients can be realized as the 
number of fix points of iterates of a set map T: X X, where the number of orbits 
of order n is divisible by n^. We figure out these very natural relations in Section 9. 

Our final table (mostly represented in the joint contribution [AESZ, Table A]) con- 
tains more than 200 cases of Calabi-Yau equations (see also van Enckevort's electronic 
database [En], which provides more advanced knowledge, like instanton numbers up 
to 20, in some cases also the number of elliptic curves, monodromy matrices etc.). 
In almost all cases, the mirror map z{q) divided by q seems to be a high power of 
a series with integer coefficients (cf. [LY]). For this we have no explanation at all, 
but we collect this experimental knowledge in another table [AESZ, Table B]. The 
last table [AESZ, Table C] presents a brief systematic guide to the main table from 
[AESZ, Table A], with entries ordered according to instanton numbers^. 

Finally, in Section 10, we find a 6th-order linear differential equation coming from 
the simultaneous approximations to ^(3) and C(5), and also a similar 6th-order equa- 
tion for the series X]^o X]fc=o (fc) • Both cases seem to admit arithmetic proper- 
ties very close to those for Calabi-Yau cases, but for the latter series we have a new 
phenomenon — a free parameter. This fact does not look apparent: we briefly indicate 
some similar well. 

Acknowledgements. First of all we want to thank Duco van Straten, who has shown 
great interest in our work and suggested the use of Hadamard products. He and 
Christian van Enckevort contributed with finding several new Calabi-Yau equations, 
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from [ES] resulted in the joint project [AESZ]. Frits Beukers gave us details of com- 
putations on [Be4] and inspiration to start this work; his comments helped us a lot to 
improve the first version of our manuscript. Helena Verrill communicated us a couple 
of new equations and Arne Meurman showed an algorithm to find the Hadamard 
product. Special gratitude is due to Noriko Yui who kindly proposed us to contribute 
the BIRS workshop volume on Calabi-Yau varieties and mirror symmetry. 
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DIFFERENTIAL EQUATIONS, MIRROR MAPS AND ZETA VALUES 



3 



1. Maximal unipotent monodromy 
Consider a linear differential equation 

+ a,_i(^)|/(^-i) + • • • + a^{z)y' + ao{z)y = 0, (1.1) 

where the prime stands for the z-derivative and ao, ai, . . . , a^-i are meromorphic 
functions in variable z. Assuming that ^ = is a regular singularity of the differential 
equation (1.1), we may write its coefficients in the form 

as^j{z) = z~^as-j{z), j = 1, . . . , s, 

where the functions as-j{z) are analytic at z = 0. The roots of the indicial equation 

A(A - 1) • • ■ (A - s + 1) + a,_i(0)A(A - 1) • • ■ (A - s + 2) + • • • + ai(0)A + ao(0) = 

determine the exponents of the differential equation (1.1) at the point z = 0. Fol- 
lowing [Mo], we will say that our differential equation (1.1) is of maximal unipotent 
monodromy (briefly, MUM) if its exponents at ^ = are all zero. The Frobenius 
method (see, e.g., [In, Section 16.1]) gives one a constructive way for writing a basis 
of solutions to (1.1) in the form 



(1.2) 



y{z, p)^Yl P^^""^' (modp^) = yo{z) + yi{z)p + ■■■+ ys-i{z)p'-\ 

n=0 

A(0,p) = 1 (modp') 
(in particular, A{0, 0) = yo{0) = 1), where by definition we set 

z'^ imodp^) = l + \ogz.p+^^.p' + ...+ ^-^^.p^-\ 

The components yo,yi, ■ ■ ■ , ys-i in (1-2) define the Frobenius basis of the MUM dif- 
ferential equation (1.1) at z = 0. 

2. 4th- ORDER LINEAR DIFFERENTIAL EQUATIONS 

Consider a linear homogeneous differential equation of order 4, 

y(^) + a3(2;)y(3) + a2{z)y" + ai{z)y' + ao{z)y = 0, (2.1) 

where ao,ai,a2;'^3 are rational functions in the variable z. Suppose that (2.1) is 
MUM and let yo,yi,y2, J/s be the Frobenius basis of solutions to (2.1). Set t = j/i/j/o- 

Proposition 1 [A12]. In the above notation, the condition 

11 3 1 

ai = -0203 - -al + 4 - -asa's - -a'^ (2.2) 
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is equivalent to any of the following two possibilities: 



ex-p(-lJ a3{z)dz) 
^(^2/yo) = and ^{ys/yo] 



d2 



Proof. It can be easily checked that the functions 

vo yo 



n = 



m 
yo 



satisfy the 4th-order hnear differential equation 

T(^) + 63^^^) + + biT' = 0, 



where 



6i = ai + 2a2^+3a3^ + 4^. 



yo 

yo yo 



yo 



yo 



62 = a2 + 3a3^ + 6^, 63 = as + 4^. 

2/0 



Since 



we have 



^(y_2 
dt^ Uo 



d_ _ /dt\ ^ d_ _ ^ _d 
di ~ [dzj dz~ ¥idz' 



2 ±(H]^]_f'^2 



dt^ dt\TiJ Ti\Ti (T^)2 



R 



(dt/dz) 



where 



R = T[T^' - T['T^ = 
With the help of (2.4) we deduce 



rpll rpll 
^1 ^2 



R' = 



Ti TI, 

rj-\lll rj-\lll 



R" = 



rpll rpll 
rplll rplll 

^1 ^2 



bsR' - b2R, 



R'" = -b^{R" + b^R' + b2R) + biR - {bsR'Y - {b2R)'. 
Therefore, the function R satisfies the 3rd-order differential equation 

R'" + C2R" + ciR' + cqR = 0, 

where 

Co = ^'2&3 -bi-\- b'2, Ci = 62 + &3 + ^3, C2 = 263. 

We now claim that the function 



1 1 /" 

^=^exp^-- j a3(a;)da;^ 
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also satisfies the linear differential equation (2.8) if and only if (2.2) holds. By differ- 
entiating we obtain 



V 2 yoj 



i^" = (M_^ + 2a3^-2^ + 6^)i?, 



4 2 yo yo Vo 

( I ^ 2\y'o 

V 2 / Vn 



R'" = ('^ - ^ - ^ + 3f4 - -ai)^ + 3a3^ -2^ 
4 2 8 V 2 /yo yo yo 



_ 9a3^ + 18^ - 24^^)^, 

yo yo yo J 

hence a necessary and sufficient condition for the function R to satisfy (2.8) is the 
required one: 

,3 



It remains to compare the first three terms (i.e., coefficients of x~^, x~'^, and x~^) in 
order to verify the equality R — R. 

We are now required to check the second equality in (2.3) provided (2.2) holds. 
Changing R to i?i and T2 to T3 in (2.6), (2.7) we obtain that the function Ri sat- 
isfies the differential equation (2.8). On the other hand, differentiating the product 
TiR = tR with the help of the differential equation (2.4) for Ti and of (2.8) for 
i?, we deduce that the function tR also satisfies (2.8). Comparing of the first three 
coefficients again yields Ri = tR, hence formulas (2.5) and 

ys\ d^Ts d fT^\ Ri 



df^yyoj dt2 dt\TiJ (dt/dz)^ 
complete the proof. □ 

If the differential equation (2.1) is MUM and assumption (2.2) holds, we define 
the mirror map z = z{q) as the inverse of q{z) = = exp(yi/|/o) and the Yukawa 
coupling as 

d2 

K{q)^No- ^{y2/yo), 

where A'^o is some integer different from 0. The interesting (Calabi-Yau) cases are 
those when z{q) G and 

Nil^q^ 

n=0 1=1 

with 



K(g) = 5:C.g- = iVo + j:^ (2.9) 



No = CoeZ and = ^ ^ Z forZ = l,2,... 

d\l 

and iJ, is the Mobius function (and, usually, Ni > 0). 
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Proposition 2. 


Under hypothesis 


(2.2) 


we 


have the 


nants: 
















yo 


y3 




yi y2 






y'o 


y's 




y'l y2 



Proof. The required identity is equivalent to VQiys/yoY = 1/1(^2/1/1)', i-e. 



t yo 



yo 



yo 



yi yo 



Since ^ = • ^, verification of the last identity is equivalent to showing that y = 0, 
where 



Y 



^fys 

dt \yo 



But 



_^2_^/^i y2 

dt\t yo 
dY d2 



d / 2/3 



dt \yo 



+ 



y2 
yo 



d ( y2 



dt \yo 



dt 



y^\-t 



d2 



^1=0 



df^yyoj dt'^Kyo, 

by Proposition 1. Therefore Y = const, and verifying near 2; = we get Y = (all 
log-terms disappear). This completes the proof. □ 



3. Wronskian formalism 

Consider an arbitrary pair y{z), y{z) of linearly independent solutions of a 4th-order 
differential equation (2.1) (not necessarily MUM). Define the Wronskian determinant 



Y 



y y 
y' y' 



Proposition 3. The function Y = Y{z) satisfies a 6th-order linear differential equa- 
tion with coefficients depending only on the coefficients 0/ (2.1) {i.e., independent of 
the choice of the pairy, y), called the exterior square o/(2.1). Moreover, this Qth-order 
differential equation reduces to a bth-order differential equation {i.e., the coefficient 
ofY^^^ is zero) if and only if condition (2.2) holds. 

Remark 1. In other words, condition (2.2) is equivalent to the exterior square having 
order 5. As pointed out by Beukers, by this equivalence properties (2.3) are the 
derived properties which are handy to calculate the Yukawa coupling. Nevertheless, 
we kept the original proof of Proposition 1 to make it independent of considerations 
in this section. 

Remark 2. The clain of Proposition 3 may be proved in Maple by the DEtools com- 
mand exterior_power. 

Proof. We have 



U4 



Y'-- 




y y 






y" y" 








y' 








1 


U5 = 


y" 







U3 



y' 



y 

y'" y'' 



y 



y 

U6 



y y 
y'" y" 
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Using the differential equation (2.1) for y, we get the system 

Y" = tt3 + tt4, tig = tt5 - aiY - a2Y' - agttg, 

tt4 = tt5, = ug + aoY - 02 - asu^, 

u'q = aoY' + aiU4 — gsUq. (3-1) 

Then 

Y'" = + ^4 = -Us - aiY - a2Y' - asUs + U5 

= 2u5-aiY -a2Y' -asiY" -u^), (3.2) 

and in the notation 

U = Y'" + a^Y" + asF' + a-,Y 
we may write (3.2) in the form 

2tt5 = U — 03164. (3-3) 

Therefore 

U' — a'2,U4^ — a^u^ = (U — 03144)' = 2u'^ = 2uq + 2aQY — 202^4 — 2031x5. 
Hence by (3.3) 

U' + ^a^U - 2aoY = 2uq + {^al - 2a2 + a'^u^ 
and after taking the derivative 

U" + -a'iU' + ^a'^U -2aQY' -2a'QY = 2uq+ ^^a|-2a2 + a3^it5 + (a3a3-2a2 + a3)it4. 

Using (3.3) and (3.1) we write the last equality as W — Vu4^, where 

W^U" + ^asU' + (as + ^4 - ^a^) V - 4aoF' - 2{o!^ + a^ai)Y 

and ^ ^ 

y = 03 — 2a2 + 770303 — 0203 + -a| + 2a\. 

The condition F = is exactly (2.2); if it holds, we have W = Vu^ — 0, which is the 
required 5th-order differential equation. If F 0, then differentiate W = Vu4 to get 

W' = V'u4 + Vu5 = V'u4 + lv{U - 031x4), 



hence 



VW = ^V'^U + (V' - ^asV^W. (3.4) 
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Equality (3.4) is the 6th-order linear differential equation for the function Y; the 
coefficient of Y^^^ in this equation isV^O. □ 

Prom now on, we assume condition (2.2) for a given 4th-order MUM differential 
equation (2.1). To construct a 5th-order MUM differential equation we modify the 
above construction by taking 



w{z) = zY = z 



y y 
y' y' 



Denote the resulting 5th-order differential equation for the function w: 

w^^'^ + hA{z)w^'^'^ + hz{z)w^^'^ + h2{z)w" + hi{z)w' + hQ{z)w = 0. (3.5) 

Remark. Clearly, condition (2.2) induces a certain condition for the coefficients of the 
differential equation (3.5). We leave the corresponding (rather boring) relation to the 
reader as a Maple exercise. 

Proposition 4. The differential equation (3.5) is MUM; its Frobenius basis is given 
by the following formulas: 



Wo = z 



yo 


yi 


, Wl = z 


yo 


2/2 


, W2 = Z 


yo 


2/3 


= z 


yi 


2/2 


















y'o 


y'l 




y'o 


y'2 




y'o 


y's 




y'l 


y'2 



W3 



2/1 2/3 

y'l y's 



W4 



2/2 2/3 
2/2 2/3 



(3.6) 



Remark. The inverse is not true in the following sense. If we take a MUM linear differ- 
ential equation of order 5 with monodromy group O5, then its 4th-order differential 
puUback is not necessarily MUM. An example to this is the 5th-order differential 
operator 

D = e^-6^z{o + l){o + l){o + l){o + l){o + l) 

(which is a natural generalization of the 4th-order operator 

D = 9'-5^z{9+l){9+i){9+i){9+l) 
corresponding to the Calabi-Yau differential equation). 

Proof. The five functions in (3.6) are linearly independent solutions to (3.5). De- 
veloping their (log z)-expansions (with the known structure of the Frobenius basis 
2/0)2/1)2/2)2/3 to (2.1)) show that they form the Frobenius basis of (3.5), hence the 
latter differential equation is MUM. □ 

Corollary (Beukers's relations). We have 

2wqW4 — 2wiWs + W2 — 0, 2w'qw'^ — 2w'iw'^ + w'^ = 0. (3-7) 



Proof. Write 



yo 


2/3 


■ z 


2/1 


2/2 










y'o 


y's 




y'l 


y'2 
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Then the first identity is trivial by expanding. In the same vein we obtain the identity 
for the functions Wj, in which the second fine of the determinant contains yls replaced 
by y'/s. Since w'j = Wj/z + Wj, we finally arrive at the second identity in (3.7). □ 

Remark. Beukers proved the above relations for a concrete example of the 5th-order 
differential equation, which we discuss in the next section. 

In the case of the 5th-order MUM differential equation (3.5) we can define analogues 
of a mirror map and a Yukawa coupling. Namely, we can define z{q) as the inverse of 
q{z) — exp{wi/wo) and 

~ (P 

K{q) ^ No ■ -^{w2/wo), t^Wi/wQ. 
Developing the last q'-series as the Lambert series, 

= ^0 + E (3-8) 

(we put instead of /"^ here), we have the following experimental observation: z{q) e 
K[[q]] and Ni e K, where K is either Z or Z[l/|)] for a certain prime p (the possibilities 
p = 3,5,7, 23 occur). 

4. Strange integrality related to C(4) 
In the work [Zul], the following 3-term polynomial recursion is given: 

(n + ifAn+i - 3(2n + l){3n^ + 3n + l){15n'^ + 15n + 4)^^ 
- 3n^(3n- l)(3n + l)A„_i = for n ^ 1. 

If we take two linearly independent solutions {An} and {-Bn} to the above recursion 
given by the initial data 

^0 = 1, Ai = 12, and Bq = 0, Bi = 13, 

then 

lim ^ ^ C(4) = - 

n— *-oo 90 

(see [Zul] for details). This recursion was previously proved in [Co] and [So] but 
without any indication of arithmetic properties of the sequences {A^} and {Bn}- 
In [Zul], it is proved that 

DnAn, D^Bn £ Z, where Dn = the least common multiple of 1, 2, . . . , n, 

and conjectured some stronger inclusions that were finally proved in [KR]. In partic- 
ular, from [KR] we have the inclusions 

An,D^BneZ, 
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and even the following explicit formula (see also [Zu4]): 



n\'^ f n\'^ f n + j\ f n + k\ f j + k 



j J \k J \ n J \ n J \ n 



where the binomial coefficients (^) are zero if 6 < or a < 6. 

A natural object related to the above recursion is the 5th-order MUM differential 
equation 

z^{27z'^ + 270z - l)w^^'^ + z^{A<dbz^ + 3375z - 10)t(;(^) 

+ ^3(1752^2 + 11502^ - 25)t(;(^) + ^^(2412^2 + 11259^ - lb)w" 

+ ^(816^2 + 2130^ - l)w' + 12^(2^ + l)w = 0. (4.1) 

Denote by wo,wi,W2,W3,W4 the Frobenius basis of (4.1) near ^ = 0, so that 

oo 

Wo = wo{z) = Anz"" = l + 12z + 804^2 ^ ggesOz^ + 12386340^^ + • • • . 

n=0 

It was pointed out by F. Beukers [Be4] that the (Zariski closure of the) monodromy 
group of the differential equation (4.1) turns out to be O5. This result (as well 
as identities (3.7) for the latter functions wjs) are consequences of the following 
statement. 

Proposition 5. The bth-order MUM differential equation (4.1) is constructed from 
the Ath-order MUM differential equation (2.1) with 

6(36^2 + 3152-1) 

as = 



02 



ai = 



ao 



^(27^2 + 270^ - 1) ' 
10530^'^ + 180387^3 + 759417^^ _ 4571^ ^ 7 
22(27^2 + 2702-1)2 ' 
962282^ + 24210092^ + 184165652^ + 493398542^ - 45005422 + 11072 - 1 

2^(2722 + 2702 - 1)3 ' 

3(87482^ + 2536922^ + 108753032^ - 376010102^ 

+136433282^ + 43213522 + 132242 + 11) 



2^(2722 + 2702- 1)4 

by the algorithm of Section 3. 

Proof. By direct computation (using Maple). □ 

Remark. The original proof by Beukers [Be5] of relations (3.7) for the differential 
equation (4.1) follows the principle 'the proof is trivial once has found the explicit 
formula'. We reproduce the arguments from [Be5] here. The symmetric square of 
the 5th order equation (4.1) has order 14, one order less than expected 15. This 
means that there should be a homogeneous quadratic relation between the solutions 
Wo, . . . ■,W4 oi (A.l). It suffices to find a relation which vanishes at 2 = of order 
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one higher than the local exponents given by the symmetric square of order 14. 
(Clearly, any non-trivial relation satisfies the 14th order equation; if we find one with 
one vanishing order higher than predicted by this 14th order equation, it must be the 
trivial solution.) If we replace the wjs by their derivatives, we get quadratic forms that 
form a one-dimensional monodromy representation. Looking at the local exponents 
we can infer that they should be rational functions with poles only at infinity. Hence 
the expressions are polynomial of degree again given by the local exponents at infinity 
of the symmetric square. 

Here are some observations related to the 5th-order differential equation (4.1) and 
its 4th-order puUback given in Proposition 5 (we denote by yo,yi, j/2, J/3 the Frobenius 
basis of the latter differential equation). 



Observation 1. We have yo e Z[l/2][[^]] but not yo e Z[[z]], while z{q), K{q) e Z[[q]] 
and for K{q) we have the Lambert expansion (2.9) with Ni e Z but Ni/Nq < 0. 

Observation 2. We have wq = Uq, where 

Uo^l + 6z + 384^2 ^ 42036^^ + 5867226;2^ + ■■■ E Z[[z]]. 

Observation 3. One has the following generalized Beukers's relations: for A; ^ 2 



where Ck,i, G Z. 

Observation 4. Put t = wi/wq, q{z) = = z + 0{z'^) and denote by z{q) its 
inverse. Then q{z) e ^[[z]] and, hence, z{q) e Zffq']]. Moreover, 



Then L„ G Z and L„ > 0. What is a enumerative meaning of the numbers (i.e., 
do they count some geometric objects)? 

Observation 6. Write 



V 



q{z) = z{l + 15z + lUbz^ + • ■ ■ )^ 

z{q) = q{l - 15g - 245g2 - 20138g2 - 2043703g^ + • • • )^- 



Observation 5. Expand wq in Lambert series: 





Then Ni E 'E and A*"; > 0. What is a enumerative meaning of the numbers Ni (i.e., 
do they count some geometric objects)? 
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Observation 7 (Kummer supercongruences) . If we write wq and K as functions 
of z, 

oo oo 

wo{z) = ^H^"' = E ^H^''' 

n=0 n=0 

then 

yl(np'") = A{np''-^) (modp^'"), C(np^) = C(np'"-^) (modp^'"), 

for primes p and (n, p) = 1; r = 1, 2, . . . . 

Recently the second author [Zu5] (see also Theorem 4 in [Zu2]), when approximat- 
ing Lis (2) (see definition (9.3) below), found two more fourth order examples related 
to simultaneous rational approximations to C(2) and C(3)- About the same time the 
first author discovered compact formulas 

^-^mnr:) 

and 

^-=[n)^[k) [n)[n)' " = 

for generating series yo{z) — ^n^" satisfying exactly the same differential equa- 

tions. These examples are cases #195 and #209 in [AESZ, Table A]. The coeffi- 
cients (4.2) is of the form similar to those in [BS2] (cf. case #26), and van Straten 
has confirmed that equation #195 indeed 'comes from geometry'. Thus we have a 
direct link between number theory and geometry. 

5. Strategy of finding new Calabi-Yau equations 

How much do we know about Calabi-Yau differential equations? We discovered a 
list of 29 interesting cases taken from different papers on the subject; these cases ini- 
tiated our systematization of Calabi-Yau equations and started the table in [AESZ]^. 
The first 14 cases in the table correspond to all known hypergeometric Calabi-Yau 
equations, while the origin of the remaining 15 known cases is due to certain vari- 
ations of multiple hypergeometric series (or multiple binomial sums). The latter 
circumstance prompted us to look on those linear MUM differential equations, whose 
analytic solutions yo{z) G ^[[2]] are single or multiple sums of certain hypergeometric 
terms. There is an algorithmic way, which is due to Gosper and Zeilberger known 
as creative telescoping, to find polynomial recursions for such sums [PWZ]. Unfortu- 
nately, the algorithm for multiple summations (being a part of the general WZ-theory 
developed by H. Wilf and D. Zeilberger [WZ]) is not implemented in Maple but solving 
auxiliary systems of linear equations is. The recursions may be easily transfcrcd to lin- 
ear differential equations satisfied by yQ{z). That was our general recipe for verifying 
known cases #15-29, and then we continued our experiments by constructing several 

^Unfortunately, we are not quite sure that the nice arithmetic properties of mirror maps and 
Yukawa coupUngs are really proved in all presented cases! In this sense, [AESZ, Table A] involves 
cases that are only expected to be 'arithmetically nice'. 
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new cases, also in combination with other hypergeometric methods — quadratic trans- 
formations (Section 6), Hadamard products (Section 7) and other transformations 
(Section 8). In several cases the algorithm of creative telescoping led to reducible 
differential equations, but then we used DEtools implemented in Maple to factorize 
corresponding differential operators. 

The way a la Apery of producing new examples of the mirror maps and Yukawa 
couplings, decribed in the previous section seems to be very occasional: we do not 
really have other examples of rational approximations to 'interesting' constants (but 
see Section 10). Nevertheless, the method of finding Calabi-Yau equations by pulling 
back some 'nice' 5th-order linear differential equations works in several other situa- 
tions as well. We stress that the 5th- to 4th-order reduction did not appear before in 
connection with problems of Calabi-Yau manifolds. 

Our table in [AESZ] extends the list of 29 known cases and contains a sporadic 
list of more than two hundred Calabi-Yau (conjecturally) differential equations. The 
table is long but we have found many more equations that we decided not to include 
(not to stretch the reader's patience too far!). In particular, we have many examples 
of type An = Efc(-1)''2""^'' (2""^) C^) ( ' ) • ^he most unfortunate thing seems to be the 
fact that we cannot see any reasonable end of our list and, therefore, we do not expect 
to be much criticized by the reader for missing a hundred or more other interesting 
cases. 

The analysis of [AESZ, Table A] shows that, besides the 14 hypergeometric func- 
tions (cases #1-14), we have a plenty of other 'simple' Calabi-Yau differential equa- 
tions of type 

9^ - ciz{ad^ + a9 + b)P(d) + C2Z^Q{e), 

where P(^) is of the form {2e + lf, (36' + l)(3^ + 2), (46' + l)(4^ + 3) or (6^ + 1)(6^ + 5), 
and Q{9) splits into linear factors of similar shapes. Here a e {2, 3, 5, 7, 8, 9, 11, 13, 17} 
and ci,C2 are usually powers of 2 and (or) 3. We also have 5th-order differential 
equations similar to the the 4th-order in 2. Their 4th-order puUbacks (derived by the 
method of Section 3) are much more complicated. Case #32 seems to be even more 
complicated: 

9^ - Zz{2d + 1)(3^2 + 3^ + 1)(15^^ + 15^ + 4) 
-3z^(^ + l)^(3^ + 2)(3^ + 4) 

(this is the differential equation that started this paper). Then there are the cases, 
when both the 4th-order differential equation and the 5th-ordcr equation (obtained 
by the wronskian construction) are very complicated. The corresponding linear dif- 
ferential operators usually appear as factors of higher order differential operators, for 
which we have no pattern at all. 

One can suspect that the equations in the list from [AESZ, Table A] are Picard- 
Fuchs equations of one-dimensional families of Calabi-Yau threefolds. In many cases, 
such as the hypergeometric ones and the Hadamard products, this should not be hard 
to see. For other cases, such as the one related to C(4), this fact will be quite non- 
trivial, but perhaps all the more interesting. This is a reason of why the equations 
are named Calabi-Yau equations. 
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Remark. We are quite surprised by the enomous amount of the 4th-order hnear dif- 
ferential equations with (expectively) nice arithmetic properties. After discovering 
equation #9, the first author for a long time thought that there were no more ex- 
amples but ^1-14. Then after finding the papers [BSl], [BS2] and examples ^^15-28 
there, it was still hard to find any new cases. The use of the Gosper-Zeilberger al- 
gorithm changed this drastically. For a while we thought that whenever was a 
binomial sum (even multiple) and it gave a 4th-order MUM differential equation, then 
we automatically got an integer mirror map z{q) and the Yukawa coupling K{q) as 
in (2.9) with integers A^;. There is, however, the following example: the generating 
series yo = Zl^o^n^'^, where 

satisfies a 4th-order MUM differential equation with polynomial coefficients of de- 
gree 8; the equation produces the non-integer mirror map z{q) (actually, z{q) e 
Z[l/2] [[g]]) but the Yukawa coupling has integers Ni. There are more examples, also 
for the cases of non-integers Ni; we indicate one in Section 7 below. 



6. Quadratic transformations and mirror maps 

A purely analytic, hypergeometric machinery — quadratic and higher-order transfor- 
mations — may be put forward for constructing new examples from the old ones. The 
method powerfully works in the case of 2nd- and 3rd-ordcr differential equations, 
when we have to substitute a suitable modular function in place of variable z; see, for 
example, [HM]. 

Our starting point is the existence of quadratic transformations for 2F1- and 3F2- 
series, like 



2^1 



a, h 
1^-a-h 



.) = (i-.)-»-.F,(5«.HI;_-; 



Az 
(1^ 



(6.1) 



due to Gaufi and 



3-^2 



1 + a — h,l + a — c 

= (l-^)-«-3F2' ' 



a. 



1 + a — h — c 
a — h, a — c 



Az 



{1-zy 



(6.2) 



due to Whipple. Substituting a = h = c = ^ \n them, we get algebraic relations 
between different mirrors that are necessarily modular maps in these cases. These 
quadratic transformations have a natural generalization to the case of higher dimen- 
sional hypergeometric series (see [Zu3], Sections 6 and 7). Similar but particular 
results are given in the following statement. 
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Proposition 6. The following quadratic transformations are available for /^F^^-series: 



4-^3 



a, 6, c, d 
1 + a — 6, 1 + a — c, 1 + a — d 



—y 

' n=0 

(c)i. (1 + a - 6 - c)n_^ 



v\ (n - zv)! (1 + a - c?)^. 



(6.3) 



-1 °° /In 

1 ^ (2a), 



4z 



E 



X 



E 

i/=0 



iy!(//-iy)!(l + a-d)j. 



(-1)^ 



(6.4) 



Proof. Writing 



4-^3 



a, b, c, d 



1 + a — 6, 1 + a — c, 1 + a — d 

(o)n {b)n (c)n 



= E 

n=0 
00 

= E 

i/=0 

00 

= E 

i/=0 



/ \n rp I 71, a + 77. 

n!(l + a-6)„(l + a-c)„,^ ' H 1 + a-d 



(-1)- 



E 



(a)^+n (&)n (c)r 



u\(l + a- d)^ ^^(n-u)\(l + a- b)n (1 + a - c), 

{a)2u{b)Ac),. 

u\ {1 + a - b)^ {1 + a - c)^ {l + a-d)/ 



-i-zr 



X 3F2 



a + 2^, b + ly, c + u 
1 + a — b + v, 1 + a — c + u 



-z . 



(6.5) 



To the latter 3F2-series we apply the quadratic transformation (6.2) and reorder 
summations to get (6.3). 

For the proof of (6.4), we proceed as in (6.5) to deduce 



3^2 



z 



a, b, c 
1 + a — 6, 1 + a — c 

°° {a)2xib)x 



E 



A! (1 + a - b)x (1 + a - c). 



A=0 



X p / a + 2A, 6 + A 



l+a-6+A 



(6.6) 



To the inner 2i^i-series we now apply the transformation (6.1) to get, as before, the 
double series for the left-hand side of (6.6). Then substituting the resulting formula 
for the 3F2-series into (6.5) we arrive at the desired identity (6.4). □ 
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Plugging ina — b — c — d— ^we obtain 

CX) X 
n— n N ^ 



n=0 



where 



A* /o \ / 



The rcsiihing series '^'^^q ^i^^ z"^ and Xl^o ""-^"^ '^^^ solutions of 4th-order 
linear differential equations. They gave two new examples (cases #30 and #31) 
of differential equations producing mirror maps and Yukawa couplings with desired 
integrality properties. It seems to be an interesting problem to find quadratic and 
high-order transformations for other cases from [AESZ, Table A]. 

Besides hypergeometric transformations, there is also an 'obvious' way to produce 
infinitely many linear MUM differential equations starting with the only one such 
example (1.1). 

Proposition 7. Let y = y(z) be a generic solution of equation (1.1) of order s with 
coefficients in C{z) and let u = u{z) be any function satisfying u' /u G C(z). Then 
the function w = y/u is a solution of the linear differential equation of order s with 
rational coefficients, and this equation depends only on (1.1) but not on its solution y . 

Proof. Substituting y — uw in (1.1) we get the equation 

(uw)'^^^ + as-iiz)(uw)'^^-'^^ + h ai{z){uwy + ao{z)uw = (6.9) 

in w unknown. Developing the derivatives by the formulas 

(uw)^"^ = Y,(^)^^^'^'w^''~'\ k = l,2,...,s, 

and dividing the left-hand side in (6.9) by tt, we obtain the sth-order linear differential 
equation for w with rational coefficients since u^^^ /u e C(2;) for any j = 0, 1, 2, . . . . □ 

By Proposition 7, taking a known example of a Calabi-Yau equation with the 
corresponding series j/o = Z^^o^^-^" ^ ^[[-^1] choosing any m e 1 + zT^z]] with 
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the property u' /u G Q(2;), we see that wq — yo/uis in ^[[z]] and w — y/u also satisfies 
a Calabi-Yau equation, with the same mirror map z{q) and Yukawa couphng K{q), 
since Wj/wq = yj/yo for j — 1,2,3. For instance, the choice u = \/l — 4pz, where 
p is an arbitrary integer, leads to an infinite family of Calabi-Yau equations for the 
sequences 




Clearly, examples of such type are not presented in Appendix A. 

We also found experimentally another (doubly) infinite classes of Calabi-Yau dif- 
ferential equations depending on positive integer parameters p and r. 

Proposition 8. Given An, one of the cases presented in [AESZ, Table A] with the 
corresponding Yukawa coupling K{q), define 




Then yo — Yl'^=o ^nZ^ satisfies a Ath-order MUM differential equation with Yukawa 
coupling K{q) = K{q^). 

Remark. It can occur that the mirror map z{q) does not lie in Z[[g]] (but the denom- 
inators of its coefficients are necessarily divisors of powers of r). 

Proof. We only give a sketch of the proof, without indicating computational details. 
Suppose that D is the linear 4th-order diff'erential operator annihilating the given 
series yQ{z) = Yl^=o^nz"' and yo,yi,y2,y3 is the Frobenius basis of the differential 
equation Dy = 0. Change the variable z i— > Z{z) = {z/{l —pz))^; then substituting 

d r d 

= z— = • Z- — 

dz 1 — pz aZ 

into the differential equation Dy = leads to a new 4th-order MUM differential 
equation that annihilates Yo{z) = yQ(Z{z)). There is no difficulty in computing the 
Frobenius basis of the new differential equation (cf. (1-2)): it is 

Yj{z) = r-^yj{Z{z)), j = 0,1,2,3. 

Therefore, T{z) = Yi{z)/Yo{z) = -t{Z{z)) and the new mirror map Z{q) is related 
to the old one z{q) in accordance with the formula 



Finally, 



d2 /y. 
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Applying now the composition 

y{z) ^ y{z) = y{Z{z)) = —^ y(( ^ ) J 

1 —pz 1 — pz \\1 — pz J J 

with the help of Proposition 7 we conclude that the differential equation for y{z) is 
of order 4 and MUM, and the corresponding Yukawa coupling is K{q^). It remains 
to verify that 

Remark. As pointed out to us by Beukers, Propositions 7 and 8 might be discarded by 
choosing a normalization for Calabi-Yau differential equations. Natural requirements 
are as follows: 

(a) the smallest local exponent at the finite singularities is zero; 

(b) the point 2; = oo is a singularity. 

However, the propositions remain useful if these conditions are abandoned. 

7. Hadamard products 

Using the algorithm of creative telescoping we found almost a hundred new Calabi- 
Yau differential equations. Then D. van Straten suggested that using Hadamard prod- 
ucts of solutions to 'nice' 2nd-order equations could give required 4th-order equations. 
Surprisingly, this was also the case for many of the examples we had found just by 
accident. 

Let 

00 00 

U=^bnZ'', V=^CnZ'^ 

n=0 71=0 

be two D-finite (i.e., satisfying a linear differential equation of finite order with poly- 
nomial coefficients) power series. Then the Hadamard product 

00 CO 

y = ^anz'^ =U*V = ^hnCnZ'^ (7.1) 
n=0 n=0 

is also D-finite (see [Sta, p. 194]). If and Dy are linear differential operators 
annihilating u and f , respectively, by Du * Dy we denote the differential operator 
annihilating (7.1). We do not know a general algorithm for computing D„ * 
but for a given pair of operators Z)„, Dy the problem is easily solved by linear algebra 
arguments. Wc have found about 30 second order MUM differential equations coming 
from binomial coefficients; Zagier's manuscript [Za] provides us with 36 second order 
examples with 'nice' arithmetic properties, although a closed binomial-sum formula is 
not known in all cases. Consider the following examples for which we were successful 
in finding the differential equations for their Hadamard products: 
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(a) ^n = 5Z(j , D = e^-z{7e^ + 7e + 2)-8z\e + lf; 

k=0 ^ ^ 

(b) = E (^) ' ^ ^) > D = d'- z{lie^ + IW + 2,)-z\e + \f; 

(c) Ar, = E (^) ' (1^) ' D = e^- z{m^ + 10^ + 3) + ^z\e + l)^; 

(d) = E (^) (t) : f ) ' ^ = - 4^(3^' + 3^ + 1) + 32.^(^ + 1)^; 



k=0 

n /f^7 \ 2 



(e) ^n = E4""'(^;!') (^^_f)' ^ = ^'-M8^' + 8^ + 3) + 256(^ + 1)2. 

Taking all possible Hadamard products of these examples (including squares) we get 
15 new differential equations. Squares will be of degree 5 and the others of degree 8. 
The products are listed in #100-107, 113-115 and 120-123. 

We may also extend the above list by the following examples (f)-(h) due to D. Za- 
gier [Za] and (i), (j) due to C. van Enckevort and D. van Straten [ES]: 

(f) ^(-I)fe3"-3M^M1|^, D = 9^-3z{3e^ + 39 + l) + 27z\9 + lf; 

(g) D = 9^- z{179^ + 179 + 6) + 72z'^(9 + if; 

(h) ^„=2ri:(-i)'(-f)(„-y:)^ 

D = 9^- 3z{189^ + 189 + 7) + 729z'^(9 + if; 
D = 6'2 - 4^(326'2 + 32^ + 13) + 4096z^{9 + if; 

a) .„=43."i:(-K-f)(:y:)^. 

D = 9^- 12z(72^2 ^ J29 + 31) + 186624z'^{9 + if. 

In case (g) , no explicit formulas for An is known. In other cases the explicit formulas 
for An were found in the following way. 
The Legendre function 



Pa{z) = 2F1 

satisfies the differential equation 



—a, a + 1 
1 



1-z 
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Therefore, the function 

is annihilated by the differential operator 

6*2 - ^(2^2 ^ 26' + + a + 1) + ^2(6' + if 

(see [Za]). We have 

a-m(i^) = i:(:)V 

\ / n=0 ^ ^ 

which gives 



where 




Putting a = —1/3, —1/4, —1/6 for a suitable choice of c we get cases (h)-(j). 
Let further 



(A) 


— 


0" 


i:» = 


6^ 


-4z{29 + lf; 


(B) 


-^n — 


(3n)! 
nl^ ' 


D = 




-3z(3^ + l)(3^ + 2); 


(C) 


-^n — 


(4n)! 


D = 




-4^(46' + l)(4^ + 3); 


n!2(2n)!' 


(D) 


-^n — 


(6n)! 


D = 


9^ 


- 12^(6^ + 1) (66' + 5) 


n!(2n)!(3n)!' 



One can prove the following formula for the Hadaniard product by the method pre- 
sented below (P and Q are quadratic polynomials and c is an arbitrary constant): 

{9^ - zP{9) - cz^{9 + if) * {9^ - zQ{9)) =9^- zP{9)Q{9) - cz^Q{9 + 1)Q{9). 

Then we get the following table of cases^ from [AESZ, Table A] corresponding to 
Hadamard products: 

^We will not write out in [AESZ, Table A] the equations in cases #133-143, since the reader can 
easily derive them himself by the formula above. We have, however, checked that one gets integers 
everywhere. 
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/AX 

(A) 


(B) 


(C) 


(D) 


(a) 


45 


15 


68 


62 


(b) 


25 


24 


51 


63 


(c) 


58 


70 


69 


64 


(d) 


36 


48 


38 


65 


(e) 


111 


110 


30 


112 


(f) 


133 


134 


135 


136 


(g) 


137 


138 


139 


140 


(h) 


141 


142 





143 



The empty set at place (C)*(h) means that the Yukawa couphng K{q) is constant. 

In our computations we apphed the following method communicated to us by 
A. Meurman. The method is, in a sense, a discrete analogue of the wronskian for- 
malism. 

Method to find Hadamard products. Assume 

= Pi(n)An+i + Qi{n)An, Bn+2 = P2{n)Bn+l + Q2{n)Bn, 
and we want to find a recursion formula for Cn = AnB^- For this, define 

Rj{n) = Pj{n + l)Pjin) + Qj{n + 1), Sj{n) = Pj{n + l)Qj{n), j = 1,2, 
Uj{n) = Pj{n + 2)Rj{n) + Qj{n + 2)P,(n), , ^ ^ ^ 
Vj{n) = Pj{n + 2)Sj{n) + Qj{n + 2)Qj{n), ^ ' ' 

and take 

T2(n) = R2S1U1V2 - RiS2U2Vi, Tz{n) = P1Q2U2V1 - P2Q1U1V2, 

T^{n) = R1S2P2Q1 - R2S1P1Q2, 

and 

Wo{n) = Q1Q2T2 + S1S2TS + V1V2T4, Wi (n) = P1P2T2 + R1R2T3 + UiU2T^. 
Then 

T4{n)Cn+4 + Ts{n)Cn+3 + T2in)Cn+2 - W^{n)Cn+i - Wo{n)Cn = 0. 
When An = Bn, the method has to be modified. Assume 

An+2 = P{n)An+i + Q{n)An, 
where P{n) and Q{n) are rational. Then 

An+3 = R{n)An+i + S{n)An, 
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where 

i?(n) = Pin + l)P(n) + Qin + 1), 5(n) = P(n + l)Q(n). 
Letting — and T{n) — R{n)Q{n) — P{n)S{n), we see that 

P(n)Q(n)Cn+3 - Rin)Sin)Cn+2 - P{n)R{n)T{n)Cn+i + Q{n)S{n)T{n)Cn = 0. 

The procedure leads usually (in cases (a)-(h)) to a 8th order differential equation 
(6th order if = Bn) that factors into a 4th order equation. Cases #167-179 were 
computed by this way. 

Actually we found four 2nd order differential equations somewhat suitable for 
Hadamard products: 

(k) A„ = 9"E(-l)'(T)(r^-1)C)' D = (,' -3420 + D-S1A0 + !)■'; 



(1) /1„ = 8"53(-1)'( D = 9^ -4.(29+1)- 64.^9 +1)^ 

D = e^- 2Az{2e + 1) - 1296^2(6/ + 1)^; 
(n) -4„=4"B-l)'tf)(,;l^:)(:). D = e^-l,AO^^f. 



Only (n) gives Calabi-Yau differential equations when we form (a) * (n), (b) * (n), . . . , 
(n) * (n) but this is just the substitution z i— > z^^'^. The Hadamard squares of (k), (1), 
(m) correspond to Calabi-Yau differential equations found from the formula 

{d^ - Az{2e + 1) - Bz^e + 1)2)*2 

= - z{2Be'^ + A^{2e + if) - Bz'^{2e + 1)(2S^2 ^ 2Be + AA^ + B) 
+ B^z^{2Be^ + 8B0^ + (4^2 + 12B)e^ + {AA^ + B)e + + 2B) 

-B^z\e + i)\ 

As a curiosity we mention that in the case (m) * (m) we have 

q = z{l + 21z + 3840^2 + • • • )12960^ 

where the series in the brackets is expected to lie in ^[[2;]]. 

Of all non-square Hadamard products only (k) * (m) gives a Calabi-Yau equation 

- 4 • 18^(26* + 1)2 - 2 • 18^^2(276)4 + 36^3 + 74^2 ^ ^qq ^ 34) 

- 32 • 18^^^(261 + 1) + 2 • 18 V (276*^ + 726»^ + 1286»2 + 72^ + 13) 
+ 4 • 18V(2^ + 1)2 - 18'2^6(^ ^ 
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To get integer instanton numbers we need to make the substitution 
then 

Ni = -48, iVs = -126, Ns = -2864, = 77958, N5 = 4942032, .... 

We can also get 5th order equations by taking the Hadamard product of (A)-(D) 
with 3rd order equations of the foUowing type: 

2 



(«)-^s(:)a-)e::r) 



k=0 

D^e^ - 2z{2e + l)(5^2 + 5^ + 2) + 64z'^{e + if; 

n \ 2 ^, \ 2 



A ^ (2ky [2n-2k\ 



A;=0 

D = e^- 8z{2e + l)(26'2 + 2^ + 1) + 256^^(6' + 1)^; 

n / X 2 / , , \ 2 



fe=0 

D = e^ - z(29 + 1)(176»^ + m + 5) + z'^{d + if; 
^"/'^ ^--^ 'n + k\{3k)l 



w ^»=i:M)'3»-(3i)("r) 



fc=0 

D^e^- z{2e + 1)(702 + 75/ + 3) + 81^2^^ ^ ;l^3. 

2 /r^r \ 2 



D = e^ - Az{2d + 1)(36'^ + 3^ + 1) + 16z^(^ + 1)^ 
(0 D = e^ - Zz{2d + 1)(3^^ + 3^ + 1) - 27z2(^ + 1)3; 
(r/) D = d^ - z{2e + 1)(11^^ + 11^ + 5) + 125^2(^ + l)^ 



(.) A„=(-64r5:(-f)(;^/^) =64"$:( 



3 ^-l/4\ Y-3/4^ ^ 



D = e^ - 8z{2e + i){8e^ + 8e + 5) + 4096z'^{e + 1)^ 

2 / /o\ 2 



J \n — k 



D = 0^- 3z{2e + l){9e^ + 96* + 5) + 729z'^{e + l)^ 

2 / r /^\ 2 



5/6 
k 
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D^e^ - 2Az{2e + 1)(18^2 ^ ^ 12,) + l^m2Az^{e + 1)= 



The last 6 examples were found in [ES] (without formulas for An). 

The general formula for the Hadamard product in this case is as follows: 

(6'2 - z{2e + i)P{e) - cz^{e + 1)^) * [e^ - zQ{9)) 
= 0^- z{2e + i)p{e)Q{e) - cz\e + i)Q{e + i)Q{e) 

and we obtain the following table of corresponding cases^: 





(A) 


(B) 


(C) 


(D) 


(a) 


39 


61 


37 


66 


W) 


40 


49 


43 


67 


(7) 


44 


53 


52 


149 


(5) 


150 


151 


152 


153 



In addition, we have another type of Hadamard products. If {A^} is given by any 
of the cases {a)-{e) with equation 

- z{2d + l)P{e) + cz^{e + l)^ 

then the generating series for the sequence {(^J^)An} is annihilated by the 4th-order 
differential equation 

- 2^(2^ + l)^P{e) + Acz^{9 + 1)2(26' + 1)(26' + 3). 

Thus we get cases #16, 35, 29, 41, and 42, respectively. 
Proposition 9. Any differential equation of type D'y = 0, where 

D' = 0^- z{20 + l){ae^ + ae + b) + cz^{0 + l)^ 

is the symmetric square of Dy — 0, where 



D = e^ - z{2ae^ + ad+^) + cz^{d + 



(7.2) 



Proof. In [A12] it is proved that the differential equation 

y(^) + S2{z)y" + si{z)y' + so{z)y = 
is a symmetric square of 

u" + pi {z)u' + pq{z)u = ^ 



(7.3) 

(7.4) 



^We computed all the other pullbacks of the Hadamard products (one of those has trivial K{q)) 
and found no new instanton numbers: they are all identical with earlier known Hadamard products 
of type (a) * (6) etc. The pullbacks are contained in the database [En]. 
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if and only if 



and then 



Here we get 



S1S2 2sl s[ s'^ S2S2 _ n 



S2 Si s\ S2 

P' = J^ ^° = T " 18 " 12- 



_ -b + cz _l-{6a + 2b)z + 7cz^ _ 3{1 - Saz + 2cz^) 

° 1 — 2az + cz'^ ' ^ 1 — 2az + cz"^ ' ^ 1 — 2az + cz"^ 

which satisfy the identity above. We also get 

l-3az + 2cz^ -2b + cz 

Pi = —r. — 7. \ Po 



z{\ - 2az + 0^2) ' 4^(1 _ 2az + cz"^) ' 

and the operator d?' /dz"^ + pi{z) d/dz + pq{z) is easily converted to (7.2). □ 
Remark. The Frobenius basis of the 3rd order differential equation (7.3) is 

2/0 = ul, yi = uqUi, 2/2 = ^ul, 
where uq, ui form the Frobenius basis of (7.4). Hence 

yo 2\uoJ 2\yoJ 2 



K{q) = ( — 1=1 = constant. 



and 

dt^ \yo 

1 /2 

The substitution 2; 1— > 42; makes also uq = y^' to have integer coefficients. Since 9 
is invariant under this substitution we get 

D^e^ - z{8ae^ + 4ae + 26) + 402^(26' + if. 

Finding such equations is equivalent to finding integer solutions to the recursion 

(n + ifun+i = 2{a ■ 2n(2n + 1) + 2b)un - 4c(2n - ifun-i 

which is similar to Zagier's computation in [Za]. 

On the other hand, we found that only 2nd-order differential operators of type 

9^ — zP{9) + cz'^iO + 1)2, where P{6) is a polynomial of degree 2, were useful for 
producing 4th-order Calabi-Yau by taking the Hadamard product of two of them. 
This, however, coincides with Zagier's list [Za] of integer solutions to 

(n + l)^^n+i - a,n{n + 1)^^ + &n^^n-i = A^n- 
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Zagier's list has 36 entries. It contains our cases (a)-(h), (A)-(C) (but not (D) 
since the coeflBcients are too large). If is a polynomial in n (that corresponds to 
8 entries in the list), then the corresponding 2nd-order differential equation factors. 
Excluding these polynomial cases, hypergeometric cases (A)-(C) and terminating 
cases, we found that certain remaining cases contain examples that we do not really 
like. For instance, the recurrence 

(n + 1) An+i — 32n{n + l)An + 256^n-i = 28^n 
admits the following integer-valued solution: 

, ^24nV/^ 1/2 yfk + l/2\_^ ,, 2k + l (2n-2kV(2k\_ 
tSn-k) \ k ) {2n-2k+ir\n-k J \kj' 

the 4th-order differential operator annihilating the series X^^g Cn) ^^z^ is as fol- 
lows: 

- 16z{2e + if {89^ + 8^ + 7) + 4096^2(2^ + 1)^(2^ + 3)^ 

and the corresponding differential equation satisfies relation (2.2); then computations 
show that z{q) ^ (and A^^ in (2.9) are not integers). This is an example promised 

in Remark of Section 5. 



8. More transformations 

Transformations of the first 14 hypergeometric cases. The success with negative 
rational numbers in binomial coefficients caused us to make the following experiment. 
Taking the hypergeometric operator in case ^9, 

consider 

Then the series yo{z) = J2'^=o -^n^^ satisfies the differential equation Dy = 0, where 

D^e^- 48z{2e + if {729^ + 729 + 41) + 2^^ • 3V(26' + 1)(26' + 3)(36' + 2)(36' + 4), 

an equation found in [ES] and mentioned as case #9* in [AESZ, Table A]. We get the 
instanton numbers 

Ni = -480, N2 = -226848, ATg = 16034720, 
N4 = 1094330202744, N5 = 4352645747040, .... 

We also note that 

2/0 = (1 + 246^ + 768132^2 + • • • )^ q = zil + 12z + 33333^^ + • • • )^^^ 
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where the series in brackets are expectively in ^[[2;]]. Interchanging 5/12 and 7/12, 
we obtain case 9**: 

2n\ ^ /^-1/12^ /^-7/12^ f-5/l2\ /-11/12 

n — k J \ n — k 



— e:)?(-f^)(-?0( 
--f:)?'-KT)t-l)' 



k 

with the corresponding differential operator 

= 9^^- A8z{2e + if (729^ + 729 + 31) + 2^^ • 3^z^{29 + if {26 + 3f- 

the couphng K{q) is the same as in #9*. 

Trying the above experiment for all equations #1-#14 we find that cases #1, #2, 
#11, #12 and #14 give integer mirror map and Yukawa coupling but there seems to 
be no way to get integer instanton numbers (by doing the substitution z — > cz). Case 
#5* has trivial K{q). We have, however, a couple of 'good' cases #3*, #4*, #4**, 
#6*, #7*, #7**, #8*, #8**, #10*, #10**, #13* and #13** (see Appendix A). This 
(together with cases #9*, #9** given above) ends the list of the 14 new equations. 

By some unusual substitions we can also consider the following additional cases. 
In case #14*, we take 

■^»--c:)?(T)(-r)rara 

and obtain the corresponding differential operator 

= 9"^ -2'^- 3'^z(29 + if (29^ + 2^ + 1) + 2^° • S'^z'^{29 + 1)(29 + 3)(39 + 2)(39 + 4). 
Then, choosing 

if h) = if (h/3)'/^) = 1 + f; 

we have _ _ _ 

Ni = -2592, N2 = -307800, = 81451104, 

N4 = 144135316512, N5 = 98667659422368, .... 

In case #2*, when 

= 9"^- 2^- bz{29 + 1)^(50^^ + 50^ + 33) 

+ 2^° • 5^^2(2^ + 1)(2^ + 3)(5^ + 4)(5^ + 6), 

we need even more drastic methods; the 'instanton numbers' land in Z[\/5]. Taking 

A-to) = irta/x/5) = i + f;f^, 

k=l ^ 
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we obtain 

Ni = -256^5, N2 = 32V5 - 35260, N3 = 1004288^5, 
N4 = 835297220, N5 = 102454248704^5, .... 

Case 7^2** corresponds to the interchange of 3/10 and 7/10, hence the change of the 
factor 50^^ + 50^ + 33 by 50^^ + 50^ + 17 in the differential operator and conjugate 
K{q). Is there a geometric interpretation of the numbers iV^? 

Differential equations inspired by Guillera's formulas. The formulas discovered 
in [Gu] caused us to try the following (again we mimic case #9). Let 

< = 3456^" V(-l)'^ f'') (V2)fe(Vl2)fe(5/12)fe(7/12)fc(ll/12)fe ^ 

k=0 ^ ^ 

Then wq = Yl'^=o -^n^^ satisfies D'w = 0, where 

D' = 9^ - 288^(20 + 1) (103680^^ + 2073606'''^ + 262944^^ ^ 159254^ + 41087) 
+ 2'^^3^z'^{9 + 1) (207360^^ + 829440^^ + 15145920^ + 1370304^ + 498143) 

- 2^^3^^z^{e + 1)(^ + 2) (2^ + 3) (2406*2 + 7206* + 793) 
+ 2^^3^^z'^(e + 1)(^ + 2){9 + 3) (360^2 + 1440^ + 1633) 

- 2^9330^5 + 1)(^ + 2)(6» + 3)(^ + 4) (2^ + 5). 

The 4th order puUback of the latter differential operator is 

D = 9^- 2^3^^(248832^^ + 414720^^ + 318528^^ + 111168^ + 14497) 

+ 2^^3^°z'^{49 + 3)(432^3 + 1116^^ + 886^ + 165) 
- 23^318^3(4^ + i)(4^ + 3) (4^ + 7) (4^ + 9) 

and, remarkably, it has lower degree than the 5th order equation (cf. case #129). We 
get the instanton numbers 

Ni = 2710944, N2 = -717640978896, = 302270555492914464, 
N4 = -171507700573958028578832, = 113303073680022744870130144224, . . . , 

and note that 

yo = (1 + 260946^ + 1405445560884^2 + • • • )^ 
q = x{l + 13295Z + 67507583411^2 + . . . )576^ j^^^^ = (1 + 56478g )^^ 

where the parentheses are supposed to have integer coefficients. 
Doing the same thing with case #3, i.e. taking 



k=0 ^ ^ 
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we get the same Yukawa couphng as for case #115 which is an Hadamard square with 



J n—k 

k 



2kV f2n-2k^ ^ 



k J \ n — k 



But the differential equations are quite different as are the solutions yo- 

We list all possible variations on the theme as cases #1-14 in [AESZ, Table A]. 

There is an even simpler way to arrive at essentially the same Yukawa coupling 
K(g) as above. Consider the 5th order hypergeometric operator 



with the solution 

(12?i)! z"" 



12/ 



- (6n)!n« ^ 

of the corresponding differential equation D'w = 0, and with the 4th order puUback 

D = e^ - 2^3^^(331776^^ + 82944^^ + 132480^ - 28224^ - 14497) 
+ 2^^3^^2(2488326'^ + 124416^^ + 25056^^ - 4176^ + 21143) 
- 2^2314^3 (-33;L77g^4 ^ 248832^^ + 60480^^ + 21312^ - 4453) 

+ 2^^?,'^'^z^d{2e + 1)(12^ + 1)(12^ + 5) 

Observe that here the puUback increased the degree from 1 to 4. The instanton 
numbers are 

Ni = -2710944, N2 = -717640301160, . . . , 
i.e. the same as above up to a sign. We also note 

yo = (1-86982^-577991455848^2 )24, q = a;(l+744U+32486988467^2^- • • )^^^ 

A third variation of the theme is the following. Consider 

<=(...r(r)|:(:)(-r)(T)(T)(-'n 

with corresponding 5th order differential operator 

D' = 9^ + 2- I2^z{2e + 1) (124416^^ + 248832^^ + 234720^^ ^ ii03046' + 21121) 
+ 112 • 12^^z^{e + l){2e + l)(26' + 3)(72^2 ^ ^^^q ^ ;l01) 
+ 16 ■ I2^^z^{2e + 1)(26' + 3)(2 + 5)(1152^2 ^ 34555, ^ 2831) 
+ 12^^ {9 + 2) (2^ + 1)(2^ + 3) (2^ + 5) (2^ + 7) 
+ 64 • 12^°^^ (2^ + 1)(2^ + 3) (2^ + 5)(2 + 7) (2^ + 9) 
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with puUback 

D = 9^ + 2^3^^(995328^^ + 497664^^ + 220608^^ - 28224^ - 35233) 
+ 2^^3^^^(51425286'^ + 5308416^^ + 2946816^^ + 97926' - 292949) 
+ 2^^3^^z^ (1866240^^ + 2985984^^ + 2011320^^ + 142920^ - 208501) 
+ 2^832°^^ (6656256^^ + 14681088^^ + 11732832^^ + 1326960^ - 1309273) 
+ 26^32^^^ (4^ - 1) (51840^3 + 160704^2 + 177908^ + 64537) 
+ 2^^33^^^(4^ - 1)(4^ + 3) (4464^2 + 13536^ + 10985) 
+ + 2) (46' - 1)(46' + 3) (46' + 7) 

+ 2^°S3^^^^(4^ - 1)(46' + 3)(46' + 7)(46' + 11). 

We have 

yo = {l + 211398;^ - 2341345644648^^ + 

and 

q = z{l- 28177^ + 245565832115^^ f^^. 

The instanton numbers are the same as above. Thus we have three fourth order 
differential equations of degree 3, 4, and 8, respectively. The instanton numbers, 
however, are invariants (possibly, with some geometric interpretation). It would be 
interesting to find the transformations between the three different solutions. 

The mirror at infinity. In [Ro] E. R0dland studied case #27 at infinity. We choose 
instead #124 from [ES]: 

D = 61^^^ - 6U(30296'^ + 5572^^ + 4677^^ + 1891^ + 305) 

+ (1215215^^ + 3428132^3 ^ 4267228^^ + 2572675^ + 611586) 

- 3^2^(39370^^ + 140178^2 + 206807^^ ^ U219W + 37332) 

+ 3^2^(5660^ + 2230^3 + 33566'^ + 22416' + 558) - 3^^z^^{e + l)^ 

where the explicit formulas for are known and one should take Nq — 61 in order 
to get integer instanton numbers (cf. (2.9)). The substitution z 1-^ 3~^z~^, y ^ zy 
results in ^ — > —^ — 1 (this works only if the highest degree term has the form (6' + l)'^) 
and we obtain the 'dual' differential operator 

D* = e^ - zibme^ + 34^^ ^ ^20'^ ^ 45^ ^ 

+ 3^2(39370^^ + 17302^2 + 22493^^ + 8369^ + 1140) 

- 3^2^(1215215^^ + 1432728^^ + 1274122^^ + 538245^ + 93222) 

+ 3^612^(30296'^ + 65446'^ + 6135^^ + 2863^^ + 548) - 3^261^2^(6' + 1)^. 

We do not know a formula for the corresponding here. 

Many differential equations with highest degree term [6 + 1)^ (in particular, all 
Hadamard products) are self-dual. There is, however, the class with 
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for r = 0, 1, 2, 3, 4, 5 (cases #22, #21, #23, #56, #71, #118, respectively). Then the 
duahty maps r 5 — r. 

Also differential equations with highest degree term (2^ + 1)^ can be reflected in oo. 
Let us take, e.g., case #55 with 

D = 99^- 12^(208^^ + 224^3 + 163^^ + 51^ + 6) 

+ 2^z^{320'^ - 9286^ - 1606^^ - 837^ - 141) 

+ 2^^z^{\m'^ + 576^3 + 467^2 ^ + 15) - 22^^(2^ + 1)^. 

Then the transformation z — > 2~^'^z~^ maps ^ to — 1/2 and gives the dual 

D* = d^ - 2^z(576^^ - 11526*^ - 724^^ - 148^ - 13) 

- 2^^22(32^^ + 992^2 - 166^2 - 570 - 6) 

+ 22^3^^(8326'^ + 7686*^ + 5566'2 + 1926' + 25) - 2'^^'i^z^(2d + \f 
Also cases #33, #99, #154 can be treated in the same way (case #154 is self-dual). 

9. SUPERCONGRUENCES AND /c-REALIZABLE SERIES 

It is time to stop and discuss the phenomenon of integer coefficients of Yukawa 
couplings number-theoretically. 

Let T: X — > X be a map of a set. Define the number of points x with period n: 

fn{T) = #{x e X : T-x = x}, 

and the number of points x with the smallest period n: 

f*{T) ^#{xeX: T^x = X, #{T'x : / G Z} = n}. 

Then /°(T) = ^/^(T) determines the number of orbits of length n. Applying the 
Mobius inversion formula to fn{T) = fdC^) we obtain 

fn{T) = lj:i^Qun (9.1) 

d\n 

Proposition 10 [PW]. Let {An}^^^ be a sequence with all ^ 0. Then, there 
exists a map T : X ^ X such that — fn{T) if and only if the numbers 

1 



n^^i^)^"^' n=l,2,.... 



are all non-negative integers. 

Proof. The '=^'-part was already done in (9.1). 

In order to prove the ' <^ '-part, we take X = {0,1,2,...} and let T : X — > X be a 
permutation consisting of Bn cycles of length n for all n = 1,2,... . Then An = fn{T) 
as required. □ 

To illustrate the above construction of T in the proof, we consider the following 
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Example. Let {^n}^i — {l, 5, 7, 17, 31, 65, . . . } be given by the recurrence = 
An + 2An-i forn = 1, 2, . . . . Then {S^j^i = {1, 2, 2, 3, 6, 9, ... } and a choice of T 
is as follows: 

T = (0) (1 2)(3 4) (5 6 7)(8 9 10) (11 12 13 14)(15 16 17 18)(19 20 21 22) ... . 
The points of order 4 are 

{11, 12,..., 22, 1,2,3,4, 0}, 

hence /4(T) = 4 • 3 + 2 • 2 + 1 = 17 = ^4- 
Assume now that for some ^ 1 we have 



^^/^(^)^deN for alln = l,2, 

d\n 



This means that the number of orbits of length n is divisible hy ^ . The sequence 
{An}n^i,2,... is said to be k-realizable if 

Bn = ^J2l^Q)^deZ for all n = 1,2,... (9.2) 

d\n 

holds. Call a power series y = X^^o ^n-^" ^ Q-ilA] k-realizable if there exists an 
integer C such that {CAn}^^i is /c-realizable. 

The most remarkable property of the examples of the Calabi-Yau differential equa- 
tions is that the Yukawa coupling (2.9) is (expected to be) 3-realizable; the pseudo- 
coupling (3.8) constructed in Section 3 seems to be 2-realizable (where in cases #15-29 
we have to replace Z in (9.2) by Z[l/p] for a suitable choice of the prime p). Another 
striking property of the first 14 (hypergeometric) cases from [AESZ, Table A] is that 
yo{z) is also 3-realizable. 

The property of being A;-realizable is related to Kummer supercongruences by the 
following result. 

Proposition 11. The following are equivalent: 

(i) {An} is k-realizable; 

(ii) Amp^ = A^pr-i (modp^^) forp prime and (m^p) = 1; 

(iii) Er=i Anz^/n"" = E^=i Uk{z^), where B n are the integers defined in (9.2) 
and 



n—l 



denotes the k-th polylogarithm. 



Proof. The proofs on pp. 3-6 in [Re] ioi k — 1 are easily generalized; note that 
Lii(x) = — log(l — x). □ 

Remark. For the first n numbers Ai, ^2, . . . , A^ the case k — 1 is well described in 
Exercise 5.2 of [St]. Then there is also an (n x n)-matrix M with integer entries such 
that Ak = Tr(M'=) for /c = 1, 2, . . . , n (see [All]). 
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10. 6th- AND HIGHER ORDER DIFFERENTIAL EQUATIONS 

We would like to start this last section by presenting two other puzzling examples 
related to 4-term polynomial recursions 

(n + l)6(41218n3 - 48459n^ + 20010n - 2871)^„+i 

+ 2(48802112n^ + 89030880n^ + 36002654n^ - 24317344n^ - 19538418n^ 
+ 1311365n^ + 3790503n3 + 46005671^ - 271701n - 60291)A^ 

- 4(2n - l)(3874492n^ - 2617900?!^ - 3144314n*5 + 2947148n^ + 647130n^ 

- 118292671^ + 115771n2 + 170716n - 44541)^^-1 

- 4(n - l)^(2n - l)(2n - 3)(41218n^ + 7519571^ + 46746n + 9898)An_2 = 

(10.1) 

for simultaneous approximations to C(3) and C(5) (see [Zu2]), and 

n{n + l)^(91n^ - 18271^ + 126n - 30)^^+1 

- n(3458n^ + 1729n^ - 2947n^ - 2295n^ + 901n^ + 1190/1^ 

+ 52n^ - 228n - 60)^„ 

- (153881n^ - 307762n^ + 185311n^ + 2960n^ - 31631n^ - 88n^ 

+ 5239n^ - 610n^ - 440n + 100)A„_i 

+ 24(n - l)^(2n - l)(6n - 5)(6n - 7)(91n^ + 91n'^ + 35n + 5)^^-2 = 

(10.2) 

for the sequence Y^^=o (2)^ (see [Pe], [SJ]). 

The differential equations for generating series of (10.1) and (10.2) are of order 9; 
nevertheless, Maple easily factorise the corresponding differential operators and we 
get the 6th-order linear differential equations 

Ad 

z^{196z - 87)3(16^^ + 752^2 _ 2368^ _ i)_ 



d 

3^^(196^ - 87)2(219522^ + 873216^^ - 27786082^ + 1235312^ + 435)^ 

+ ^2(196^ - 87)(858981762^ + 2768881024^^ - 88281697562^ 

M 

+ 7144975624^2 - 1768825884^ - 491985) — 

dz^ 

+ 2^2(209321100802^ + 5086094003202^ - 16135727761442^ + 17390406950002^ 

- 8242645169042^ + 1481959336322 + 29632635)-^ 

d2'^ 

+ 2(355092917762^ + 5797121917442^ - 15303515853922^ + 15009935198242^ 

- 7316582974562^ + 1732520938862 + 20413593)-^ 

d2^ 
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+ (7138000128;2^ + 55844570880^^ - 123313425872^^ + 96429989856;23 

- 16021623504^2 + 17983065996^ + 658503)-^ 

dz 

- 6(15059072^^ - 4148928^^ + 9924264^^ ^ 214891044^^ 

+ 106071966^ + 4609521) = (10.3) 



and 



z^{z - 1)(64^ - 1)(27^ + 1)(75^3 + 1420^^ + 561z + 9)- 



dz^ 

+ ^3(2786400;?^ + 52926750^^ - 22883417^^ - 19427551^^ 

- 654306^2 + 1308^ + 126) -p^ 

^dz^ 

+ 5^2(3822480^^ + 77258112^^ - 18886036^^ - 23855750^^ 

M 

- 758207^2 _ 12752 + 99) — 

dz^ 

+ 5^(9906840^^ + 214908768^^ - 10697346^^ - 53272456^^ 

d3 



1605389^2 - 7047^ + 117)- 



dz^ 

+ 2(223704002^ + 528293250^^ + 82763885^^ - 96071120^^ 
- 2836582^2 _ 15165^ + 72)- 



dz"^ 

+ 6(16845002^ + 445177002^ + 171024752^ 

- 46873352^ - 1589672 - 540)-^ 

d2 

+ 12(157502^ + 5031752^ + 3272052^ - 18452 - 1044) = 0, (10.4) 

respectively. (The corresponding recursions will be 7-term but we do not use this 
fact.) 

The differential equation (10.3) is MUM, so taking the Probenius basis yo, yi, y2, 
. . . , j/5 we may compute the inverse of q{z) = exp(j/i/j/o), 

z{q) + 230^2 _ 26827^^ + 24147708^^ - 23334210874^^ + 26542920855790g^ 
- 33500728089853156g^ + 45492345805504886104g^ + O(g^), 
and pseudo-coupling 

where 

Ni = -320, N2 = 118264, N3 = -84117876, N4 = 80349364184, 
TVs = -90632838175404, Nq = 113783008482427048, 
N7 = -153937885949108788148, Ng = 220092295805975113694144, . . . 
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are expected to be all integers. 

The case of the (not MUM) differential equation (10.4) is more delicate: the ex- 
ponents at the point z — are 0, 0, 0, 0, 0, 1, so that we get two linearly independent 
analytic solutions at the point. The choice of the first three components of the basis 
ai z — are as follows: 

yo = l + 2z+ 66z^ + U60z^ + 54850^^ + 2031252^^ + 86874564^^ 
+ 3848298792^^ + 180295263810^^ + 0{z^) 

oni 2 r..A 3 392153 4 37695381 5 
yi = yo log z + 6z + 201z^ + 5114^"^ H — z^ H z^ 

1636192577 « 515198795302 7 19496105427567 . q, 

+ z^ + z^ + z^ + 0(zn, 

5 35 28 ^ " 

log^z , 45 2 7085 o 1750325 , 293631899 . 

= -yo^ + y^ logz + p^z+ -z^ + ^ ___,4 ^ ^^.^ 

4039916881 . 4451078944741 7 663063626523441 

+ ^7^:^ -2 + 77z;z:^ + 



200 4900 15680 

To 



, 3138305727 
+ 108^2 + 3690^3 + 164379^^ + 6851805^^ + z^ 



73262810391 7 99760681900977 g, 

+ z^ + z^ + 0(zn 

5 140 ^ ' 

(for the series j/2 we have an additional 'freedom' /3 e Q). Then we get the expansion 
z{q) = 9-69^-13593_380^4_24ggQ^5_ggg3gg^6_2gi5330257_gQ^ggg^Q4^8^Q^^9^ 

and the coefficients in the Lambert expansion (10.5) 

iVi = iVa = 3(/3 + 30), iVg = 63(/? + 20), N4 = 1357/3 + 43200, 
N5 = 5(9139/3 + 235152), Nq = 3(479743/3 + 14248200), 
N7 = 126(419883/3 + 11883130), iVg = 3(653810477/3 + 19370535600), . . . 

are expected to be in Z/3 + 

We have several further examples of 6th-order differential equations with the same 
collection of the exponents at the point 2; = as for equation (10.4). In all such exam- 
ples we get the similar phenomenon of admitting a free parameter /3 for the coefficients 
in the pseudo-coupling (10.5). Here we only mention formulas for {An}n=o,i,... in the 
examples, where the series yo = X^^^ Anz"^ form analytic solutions in ^[[2]] to the 
corresponding differential equations: 

^„=(5ff (")*, .„=(i#5:('!)(")V* 

Similar phenomena seem to be happened in more general situations. For example, 
the series 

n / X 4 



-w-E-^E(: 



n=0 fe=0 
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gives an analytic solution to the 8th-order linear differential equation Dy = 0, where 

D = e'^ie - 1)2 - 16ze^{2e + if {46 + 1)(4^ + 3){36^ + 3^ + 1) 

- 256^2^20 ^ i^^2e + 3){4e + l){4e + 3)2(4^ + 5^(49 + 7); 

a suitable choice of the solution yi{z) leads to the inversion z{q) of the series q{z) = 
exp(|/i/?/o) satisfying (conjecturally!) z{q) G Z[[g]]. Finally, one can take a free pa- 
rameter /3 for the solution y2{z) giving the pseudo-coupling (10.5) with the coefficients 

Ni = /3, N2 = 199/3 + 308464, N3 = 17(18671/3 + 18181888), 
iV4 = 2(233056091/3 + 290786122832), N5 = 3(311572981529/3 + 381333083404544), 
Nq = 2020548109265033/3 + 2585332849682835728, . . . 

that are expected to be in Z/3 + Z. Another interesting example is the series 



n=0 k=0 



n / X 7 

n 



which also satisfies a (linearly irreducible) 8th-order differential equation Dy = with 
exponents 0, 0, 0, 0, 0, 0, 1, 1 about z = 0. 
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